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INTRODUCTION 


The equation 

y" + ay' - ay + ty 3 = 0 

has been proposed by H. Lashinsky 1 as an analog to the Van der.Pol 
equation 2 to describe nonlinear saturation of aperiodic instabilities 
such as the Benard instability in fluids 3 and related plasma instabilities 4 
In the Van der Pol equation, the saturation is described by a 
nonlinear dissipation term of the form (- o' + |3y + Y y 2 )y’. The 
instability corresponds to the negative sign of the coefficient a: 

This negative dissipation feeds energy into the system. The nonlinear 
terms saturate the instability. 

In the Lashinsky equation, on the other hand, the saturation is 
described by a nonlinear restoring force of the form (- a + by 2 )y. 

The instability itself corresponds to the negative sign of the coefficient 
a in the restoring force. The dissipation term consumes energy in this 
case . 

OSCILLATORY SOLUTIONS 

First we wish to determine whether equation (l) has periodic or 
oscillatory solutions for a > 0, o' > 0, b > 0. The equation (l) is 
a special form of the more general Lienard equation 

y" + f(y) • y’ + g(y) = 0 
which has been investigated extensively 5 ’ 6 . 
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The Lienard Theorem 6 * 7 * 8 gives sufficient hut not necessary conditions 
for the existence of a strictly periodic solution y(t). The conditions 
to he satisfied are: 

1 . f(y) = a, g(y) = - ay + hy 3 analytic, which is satisfied for (l). 

2 . gCy) an odd function, g( 0 ) = 0 and yg(y) ^ 0 , so that a restoring 
force exists. This condition is also satisfied, depending on the 
values of the coefficients a and h. 

3 - f(y) an even function and f( 0 ) < 0 (negative damping at the origin, 
so that there is an unstable singular point at the origin). For 
( 1 ) f( 0 ) is positive, however, and the damping term is not a source 
of energy. 

,-y fir a 4 

= Jo f(y)dy = °W> = Jo^ (y)dy = - g ^ + b 4 go to co as 

y aa, which insures periodicity. This condition is satisfied. 

5- F(0) = 0 has a unique root at y , y^ > 0 and is a monotonically 

increasing function for y = y Q . Since we have y Q = 0, the condition 

y Q > 0 is. not satisfied. 

Since all the conditions for periodicity are not satisfied, the 
Lashinsky equation cannot possess a periodic solution. Also the Bendixson 
Theorem 5 , the Frommer Theorem 8 , the Levinson- Smith Theorem 5 * 6 show that 
(l) has no periodic solutions in the sense 

y (t + t Q ) = y (t ) ( 3 ) 

This result, however, does not preclude the existence of oscillatory 
solutions, i.e. solutions with several maxima and minima. Apparently 
these oscillatory solutions are damped. 


In order to investigate these solutions, we use the phase plane 
method 8 . We write (l) in the form of two first order equations 

y ' = u - P (u , y) 
u ' = - Ck'U + ay - by 3 = Q (u,y) 

Critical (equilibrium) points are then given by y ' =0, u' = 0; 
these points are i.e. 



so that y" = 0 at these points. The characteristic equation that 
determines the nature of these critical points is then 

\ 2 - (P + Q )\ + P Q - P Q =0 
y W y U U y 

From (4) we have P y = 0, P y = 1, ^ = a - 3by 2 , = - a, so that 

(6) becomes 

X 2 + aX - a + yoy 2 - 0 


which gives for the point y = 0, u = 0 


zb , / *r — + a 


a 

2 V if- 


so that the origin of the phase plane is a saddle point (two real roots X 
with opposite signs), so that this solution is unstable and y ro . 
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/ EL 

u -* oo for t - 00 . The critical points y = ±^j — , u = 0 give 

*i,2 = " f ± - 2a (9) 

so that these points are either a stable focus (complex ^x }2 with negative 
real parts) for a 2 < 8a, indicating a damped oscillation, or a stable 
mode (both X are real and negative) for a 2 > 8a, indicating that no 
oscillation occurs and y -* const or zero, u -* 0 for t -* «>. 

In the case given by Lashinsky 1 , o' s=» 10, a « 1 and we have no 
oscillating solution, see section 3- In this section we investigate 
oscillatory solutions, assuming a » 0, 1, a « 1. It is of interest that 
the value of b does not enter directly into this discussion which is 
based on a linearized stability analysis. This feature is noteworthy, 
since both terms ay' and - ay are nonoscillatory and the only term that 
produces oscillations is by 3 . The coefficient b enters the discussion 
only through ( 5 ). 

In order to eliminate b completely, we make the transformation 

y(x) = a 1 ^ 2 b 1 ^ 2 u (t ) (10) 


which gives 


u " + au ’ - au 4- au 3 = 0 


( 11 ) 


In the oscillatory case a 2 < 8a and the smaller a the weaker the 
dissipation. Neglecting au ' as a first approximation, by multiplying (ll) 
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by u' and integrating we find 


u ' 2 au 2 au 4 

2 " 2 + k 


( 12 ) 


where E is an integration constant (the energy). In this conservative 
system no energy is lost and energy oscillates between kinetic and 
potential energy. If a small dissipative term cm ' is now introduced, 
then E decreases slowly and E = const must be replaced by E(t). This 
procedure can be used to establish a perturbation method based on the 
averaging process 9 . 

From (12) we have 


t (u ) + cp = j 


du 




2E + au 2 - 


(13) 


where cp is another integration constant. The integral is an elliptic 
integral and may be solved by standard methods 10 . 

The inverse function u(t) - and also y(t) - may be expressed by 
Jacobi elliptic functions whose modulus k depends on the energy constant 
E lljl2 . : 3f E varies slowly because of the dissipative term cm', then 
k -* k(t) and we can use the averaging method we have described elsewhere 9 . 


6 


WOW- OSCILLATORY SOLUTIOWS 

This is the case discussed by H. Lashinsky 1 . Ebr a « 10, a « 1 
the non-oscillatory condition o' 2 > 8a is satisfied. We start from (ll) using 
the; transformation art = ? , so that 

u + u - £u + £u 3 =0 (14 ) 

where u = du/d? and & = a/ar 2 = 10 2 . Wow the perturbation problem is no 
longer singular (as it was for o' 2 < 8a) and we may set up a Poincare 
expansion of successive approximations 

u(§ ) = u q (§ ) + £ m(§ ) + £ 2 u 2 (§) + •• . (15) 

Substituting into (l4) we obtain 


^0+ Uo= 0 

(16) 

iii + Ui = u - u 3 

( 17 ) 

00 


U 2 + u 2 = U! - 3 u 0 2 Ul 

( 18 ) 


o ) 

neglecting higher orders than £ . Integration yields 
u Q (g ) = Ae“ ? + B 


( 19 ) 
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ui(S ) = - Af e 


- A 3 6 e 


.39 


3A g B 


’ 2S + 3AB 2 ie _i 


(-B+B 3 )§ + Ce s + D 


( 20 ) 


For ^ -* the terms containing B 3 | tend to infinity, so that the 
expansion (15 ) does not converge for a finite number of terms. These 
terms are typical secular terms. In the oscillatory solution the secular 
terms are avoided by the averaging process 9 . In the present case the 
secular terms are eliminated by appropriate adjustment of the arbitrary- 
integration constants (Lindstedt method 13 ). To do this we put B = 0, 

D = 0. The terms §e”^ etc. are not secular. For u 2 (?) we have 


*2. 


> -TS 


A e 


-33 


\ A 3 §e 


.3§ 


— A 5 e 
40 A 6 


_5i> 


+ 


§ 2 e 



A 2 C -3? + Ee“^ + F 

— 6 


( 21 ) 


For § -* ® we now have from (15 ) and ( 19 ) - (2l) 


u (°° ) = £ 2 F, u 1 (°° ) = 0, u'^ 00 ) = 0 


( 22 ) 


so that from (14) F = £ 2 follows. 

The following two initial conditions are of physical interest: 
u(0) = d, u(0) = f 


where 

a) f = 0, d given; b) f given, d = 0. 

We then have two equations for the three remaining constants A, C, E so 
that one constant may be chosen arbitrarily e.g. to give a better fit 
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to data obtained from a numerical integration. However, in order to adjust 
a solution of the type ( 15 ) to the solution obtained from a numerical 
integration, higher powers of £ than the second must be taken into 
account. Thus a straightforward numerical integration is simpler than the 
expansion ( 15 ). The expansion, however, gives more physical insight into 
the phenomena. The energy fed into the system at t = 0 is used to excite 

several modes - the higher the initial energy, the greater the number of 

modes that are excited. The higher the order of the mode, the more rapidly 

it decays, so that after a certain time only the asymptotic value 1 remains. 

“2 

The result of the numerical integration for £ = 10 , f = 0, 

d = 0.2 are shown in Fig. 1; the values are given in Table 1. For £ = 10 2 , 

f = 0.2, d = 0 we find the values in Fig. 2 and Table 2. 
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Table 1 


Table 2 


X 

y 

X 

y 

0 

0.20000 

0 

0 

O.09999 

0.20009 

0.04999 

0.009754 

0.49999 

0.20020 

0.09999 

0.019033 

0.99999 

0.20070 

0-99999 

0.126628 

4.99986 

0.20776 

4.99986 

0.20461 

9.99969 

0.21783 

9.99969 (l) 

0.21584 

19.9973 

0.23928 

19-9973 

0.23713 

39.9863 

0.28779 

39.863 

0.285267 

59-9757 

0.34410 

59-9753 

0.34120 

79-0148 

0.40493 

79.0148 

0.40169 

100.003 

0.47911 

100.003 

0.47556 

120.992 . 

0.55836 

120.992 

0.55466 

126.988 (a) 

0.58143 

127.988 (a) 

0.58156 

139.981 

0.63125 

139.981 

0.62758 

160.020 

0.70539 

160.020 

0.70149 

180.009 

0.77259 

180.009 

0.76956 

200.998 

0.83272 

2OO.998 

0.83022 

219.987 

0.87670 

219.987 

0.87471 

250.021 

0.92699 

250.021 

0.92570 

280.004 

0.95840 

280.004 

0.95763 

309.989 

0.97684 

309.989 

0.97639 

343.984 

0.98959 

549.984 

0.98939 

375.912 

0.993722 

375-012 

0.99360 

400.000 

0.99621 

400.000 

0.99614 

(a) Here we have 

(b ) Hear here we 

a maximum of y ' . 
have a minimum of y ' . 








Figure 2. 


